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We present a derivation of the Redfield formalism for treating the dissipative dynamics of a time- 
dependent quantum system coupled to a classical environment. We compare such a formalism with 
the master equation approach where the environments are treated quantum mechanically. Focusing 
on a time-dependent spin-1/2 system we demonstrate the equivalence between both approaches by 
showing that they lead to the same Bloch equations and, as a consequence, to the same characteristic 
times Ti and T2 (associated with the longitudinal and transverse relaxations, respectively). These 
characteristic times are shown to be related to the operator-sum representation and the equivalent 
phenomenological-operator approach. Finally, we present a protocol to circumvent the decoherence 
processes due to the loss of energy (and thus, associated with Ti). To this end, we simply associate 
the time-dependence of the quantum system to an easily achieved modulated frequency. A possible 
implementation of the protocol is also proposed in the context of nuclear magnetic resonance. 
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I. INTRODUCTION 

The rapid development of quantum information sci- 
ence has brought together several areas of theoretical 
and experimental physics pp. Much effort has been con- 
centrated in the search for solutions to sensitive prob- 
lems that prevent the efficient realization of quantum in- 
formation processing [2]. We first mention the system- 
environment coupling which induces the decoherence of 
quantum states pj], apart from other barriers such as 
scalability [J] and optimal control of individual systems 
[5]. These challenges motivate both fundamental physical 
phenomena and outstanding technological issues such as 
individually addressing quantum systems, separated by 
only few nm, with small errors [B] . 

Potential platforms for the implementation of quantum 
logic operations appeared in many fields such as con- 
densed matter, quantum optics, and atomic physics pQ. 
However, the problems mentioned above are faced by all 
the different communities when employing their particu- 
lar techniques. In the particular case of the dissipation 
and decoherence phenomena — in which we focus in the 
present work — the Redfield formalism [7J |S] and the mas- 
ter equation [S] have been the most applied approaches 
to address the environment effects on the proposed pro- 
tocols for quantum information processing. Whereas the 
semiclassical Redfield formalism relies on a classical noise 
source, a quantum environment is assumed in the master 
equation approach. In this article, considering the gen- 
eral case of a time-dependent system, we discuss general 
similarities and differences between both approaches and 
show that they are equivalent, in the sense that they lead 
to the same phenomenological Bloch equations [10 . Con- 
sequently, both of them result in the same characteristic 
relaxation times T\ and Ti associated with the longitudi- 



nal and transverse relaxations, respectively 0[H]. From 
this identification we show how these characteristic times 
are related to the operator-sum representation [2] and the 
phenomenological-operator approach 

The Redfield formalism was intended to offer a micro- 
scopic description of the relaxation phenomenon, thus 
providing a deeper understanding of the parameters T\ 
and T2. Whereas the classical noise source employed 
in the Redfield theory suffices to derive both relax- 
ation times, two distinct quantum environments must 
be adopted to derive these time scales from the mas- 
ter equation formalism. On this regard, an amplitude 
and a phase damping environment are assumed to de- 
fine the longitudinal and the transverse relaxation times, 
respectively. These quantum environments represent an 
energy-draining and a phase-shuffle channel by which the 
system loses excitations and phase relations. 

After presenting a detailed derivation of the Redfield 
theory and comparing the derived characteristics times 
with those obtained from the master equation, we fi- 
nally apply these equivalent formulations to the prob- 
lem of state protection. We note that several distinct 
techniques have been proposed to control the effects of 
decoherence on quantum states, aiming to enlarge the fi- 
delity of quantum information protocols. Among others, 
we mention the quantum-error correction codes |12j . en- 
vironments engineering [13| . decoherence-free subspaces 
[14] . and dynamical decoupling [15]. We finally men- 
tion that in a previous work [16] . addressing the energy 
draining and decoherence of a harmonic oscillator, it was 
demonstrated that the inevitable action of the environ- 
ment can be substantially weakened when considering 
appropriate non-stationary quantum systems. Reason- 
ing by analogy with the technique presented in Ref. [16] , 
we show how to enlarge the longitudinal relaxation time 
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associated with the amplitude-damping channel focusing 
a spin-1/2 system. The ideas presented here for deco- 
herence control can be easily implemented in the nuclear 
magnetic resonance (NMR) context. 

This article is organized as follows: In Sec. II we 
present the derivation of the Redfield equation for the 
general scenario of a time-dependent system. In Sec. 
Ill we apply the master equation approach to the same 
case. In Sec. IV we show the equivalence between the 
Redfield and the master equation formalisms by deriv- 
ing the Bloch equations from both approaches. Focusing 
on a time-dependent spin-1/2 system, in Sees. V and 
VI we present the operator-sum representation and the 
phenomenological operator approaches and their relation 
with the previous techniques. As an application of the 
theory, in Sec. VII we address the state protection of 
a non-stationary spin-1/2 system. Finally, Sec. VIII is 
dedicated to our final remarks where we discuss the gen- 
eralization of the methods presented in this article for 
larger systems. In general, we will adopt the language of 
the NMR quantum information processing [17] , although 
the theory presented here is valid for several other plat- 
forms, as quantum dots [18] . superconducting artificial 
atoms [Tjj5J, etc. Throughout the article we will use nat- 
ural units such that h= ks = 1 • 



II. REDFIELD FORMALISM FOR A 
TIME-DEPENDENT SPIN SYSTEM 

Considering the interaction of a time-dependent spin 
system, described by the Hamiltonian Hs(t), with a spin 
lattice modelling the environment and represented by the 
Hamiltonian Hl, the total density operator p sch (t) in the 
Schrodinger picture, evolves as 



dp sch (t) 
dt 



-i[H s (t) + H L + H SL (t),p sch {t)], (1) 



HsL{t) being the time-dependent spin-lattice interac- 
tion. In the interaction picture, defined by the unitary 
transformation U(t) — exp {— i (7is(t) + Hrf)}, where 
Hs{i) — J dt'Hs(f), we simplify the above evolution 
equation to the form 



these pulses, the prepared state p sch (0) of the system, de- 
scribed by the diagonal Hamiltonian Hs{t), evolves only 
under the action of the environment. 



By its turn, the density operator in the rotating 
frame is given by pit) = e l ' H s{t)'pSch^ e -iHs{t) wnere 
p sch (t) = e lHLt p sch {t)e~ lHLt . By assuming a weak 
system-environment coupling and getting rid of the de- 
grees of freedom of the spin lattice, we solve Eq. ([2| up 
to second order of perturbation theory, obtaining 



da(t) 
dt 



-Tr 7 



dt' [V SL (t),[VsL(t'),p(t')}} 



We considered that the interaction Vsx(i) is a stochas- 
tic operator with null mean value [8] , which results in 
a zero first order term. Next, let us use the Markov 
approximation p(t') —> p(t) ~ a{t) ® Pl(0), <r(t) be- 
ing the reduced density operator of the spin system, 
a(t) — Ti~Lp(t), and ph the reduced density operator 
of the environmental spin lattice. This approximation 
means that the state of the lattice is not affected by 
the interaction with the system. In other words, it 
means that the lattice presents a sufficiently large heat 
capacity in order to remain in the thermal equilibrium 
state p L (0) = e-P HL /Tr [e~ pHt -], with /3 = l/T, T be- 
ing the environment temperature. Finally, inspired in 
NMR systems [17], we are going to apply the high- 
temperature approximation, which takes into considera- 
tion systems were the energy gap between the spin lev- 
els, tuv (where u) is the characteristic transition frequency 
among levels), is much smaller than the thermal energy, 
fcflT, of the system, i.e. hio/hsT « 1. In this sense, 
the density operator of the system can be written as 



a(t) 



-PH 8 (t) 



/Z ~ 1 - /3H s (t). We stress that there 



is a crucial difference between high and infinite tempera- 
ture limits; differently from the latter case, in the former 
there is still a population difference between the spin lev- 
els which accounts for the reminiscent equilibrium mag- 
netization. Thus, applying the high-temperature approx- 
imation we obtain a(t) £g> Pl(0) — cr{t) — fiH^ and, con- 
sequently 



dp(t) 
dt 



-i [V SL (t),p(t)}, 



(2) 



where V SL {t) = e ms ^nsL(t)e~ ms ^ , with H SL (t) = 
e lH Lt]j SL (t^ e -iH L t yy e nave assum ed the condition 

[Hs(t), Hs(t')) = which is always fulfilled whenever the 
free Hamiltonian of the system can be written as a diag- 
onal time-independent operator, with a time-dependent 
coefficient Hg (t) = f (t) Os- Moreover, we observe that 
in NMR relaxation experiments all the required pulses 
to perform the necessary rotations are applied either at 
the beginning, to prepare the initial state p (0), or at 
the end of the experiment, to implement the tomogra- 
phy of the evolved state. Between the applications of 



^ = Tr L f* dt' [V SL (t), {Vs L (t'),(3H L - a(t)}] . (3) 



From the Heisenberg equation of motion for Vsl(£), we 
obtain 



[V SL {t),a(t)-PH L }=[Vs L {t) >( T(t)+pHs{t)} 

dt 



■i/3W(t) 



fu ,dH SL (t) 



dt 



U(t), 
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and, consequently, the evolution equation reads 
da(t) 



dt 



= iPTr L {[V SL (t),V SL (0)}} 



Tr L / dt' [V SL (t),[V SL (t'),a(t) + /3H s (t')}} 

r v SL (t),uHt') dH ^lu(t') 



i/3Tr L [ dt' 
Jo 



(4) 



By rewriting the spin-lattice interaction as Vsi,(t) oc 
X(t)0, where X(t) models the lattice stochastic fluctu- 
ation and O stems for an operator acting on the spin 
system space, we verify straightforwardly that the first 
term of the r.h.s. of Eq. Q is null, in accordance 
with the assumption (A(t)) = 0. Moreover, with the 
above definition for the spin-lattice interaction, we ver- 
ify that HsL{t) = HsL(t) and, consequently, Vsx(t) = 
e *Hs(t) HsL fy e -iHs(t) _ Integrating by parts the third 
term of the r.h.s. of Eq. Q and considering, as usual, 
that the time oscillations of the operator W (t')OU(t') 
is much faster than that of Vgz,(i), we apply the ro- 
tating wave approximation to conclude that this term 
is also null. The fact that this is indeed the case, 
can be seen as follows: the operator VsL,{t) oscillates 
with the spin-lattice coupling frequency, while the oper- 
ator W (t')OU(t') oscillates with the bare spin frequency 
which (in the assumed system-environment weak cou- 
pling regime) is much higher than the interaction fre- 
quency. Putting all this together, we finally obtain the 
simplified equation of motion for the spin system 

^ = Tr L J* dt' [V SL (t), [V SL {t'),pH s {t') - a(t)]} , 

which, in accordance with the high-temperature approx- 
imation, where a eq (0) ~ 1 — /3Hs(0) and /3Hs(tf) — 
1 — cr eq (0), becomes 



dE(t) 
dt 



dt'[VsL(t),[Vs L (t')Mt)}}, (5) 



where we have defined the operator E(t) = cr(f) — a eq (0) 
and substituted the trace over the lattice degrees of free- 
dom by the ensemble average over stochastic realizations, 
represented by the over bar. We have thus obtained the 
Redfield equation for a time-dependent spin system and 
we note, in spite of the c-number character of the en- 
vironment degrees of freedom, its resemblance with the 
master equation to be presented below. 

It is important to stress that the high-temperature ap- 
proximation, allowing us to define the latter operator 
E(t), indicates a relaxation to the highly mixed thermal 
Gibbs state. However, we mention that in the whole cal- 
culation to obtain the Redfield equation ^ it is not nec- 
essary to impose such approximation. It was only done 
because it is characterisitic of NMR systems, on which 
we focus in the present work. 

Towards the definition of the lattice spectral den- 
sity we next introduce, through the eigenvalue equation 



Hs(t)\k) = €k(t)\k), the spin basis {\k}}. Taking the 
matrix element kk' of Eq. (pi) and back to Schrodinger 
picture where Z sch (t) = e -** s W£(*)e ms tf>, we obtain 
the Redfield equations for the evolution of the density 
matrix elements 



dt 



= -i(k\ [H s ,X sch ] \k') 

+ e" i(n **' +n "'- ) i2fcn,n'*'E^, (6) 



where we have used the short-hand notation Cl kn (t) = 
E k {t) - E n (t), E k (t) = f*dt'e k (t'), E$(t) = 
(n\Yi sch (t)\n') , and denned the relaxation matrix ele- 
ments 

■^kn,n'h / (J*) — Jkn^n'k' (^s^n'fe')^ ^ ^ 

J n 1 k' .kn(t") ^fcn)^ n k ) 



^k'n' ^ ^ Jkj.jni^-) ^jn) 



&kn ^ Jjk',n'j i^"> ^n' j)& ^ k ' j (7) 



with the environment spectral densities given by 

A„,n'fc'0,fW)= / dt'G kn , n , k ,{t,t'y n ~'« l{ t'\ (8a) 
Jo 



G kn , n > k ,(t, = (M H SL (t) \n) (n'\ H SL (t') \k>). (8b) 

To simplify the notation we have omitted the explicit 
time dependence of all functions in Eq. (JsJ) . 

For the particular case of a time-independent system, 
we obtain E m (t) = e m t, Q kn (t) = (e k - e„) t = u) kn t and 
the above Redfield equations reduce to the well-know text 
book result E] 



dt 



-i{k\ [H s ^ sch (t)] \k>) 

+ E e- i(w "'-^' Jt ^,n' fe '(*)^(i). 



with 



tt'kn,n'k' J kn,n' k' y^n' k' * J n' k' ,kn\uJkn)& 

&k'n f ^ Jkj (^j'ti)^' j 
3 

— S kn E Jjk',n'j{Unj)e lUjklt i 



and 



Jkn.nn'i^nn') — / dt'G kn , n n' (t ) exp {lUJ ml d'} . 

Jo 

We observe that, although we have focused on a spin 
system, the equations obtained here are completely gen- 
eral, being valid for whichever the Hamiltonian Hs(t), 
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provided that the three following conditions are met: i) 
[Hs(t),Hs(t')] = 0, ii) system-environment weak cou- 
pling regime (Markovian environment), and iii) high- 
temperature approximation (specifically, to derive Eqs. 
Q, ([7]), and ([8])). The restrictions and the validity of 
these approximations will be discussed in the conclusions 
of the article. Let us now turn to the master equation 
approach. 



III. THE MASTER EQUATION APPROACH 

In this section, in contrast to the semiclassical ap- 
proach of the Redfield formalism, we derive the master 
equation governing the dynamics of the dissipative time- 
dependent spin system where the environment is assumed 
to be modelled within the quantum formalism. We start 
from Eq. pL such that 



da(t) 
dt 



-Tr, 



dt'[VsL(t),[VsL(t'),cr(t)]]. (9) 



Now, instead of assuming a classical environment leading 
to the above defined spin-lattice interaction as VsL(t) oc 
X(t)0, we consider two distinct quantum environments, 
to be defined below as the amplitude- and the phase- 
damping channels, each one being modelled by an infinite 
collection of decoupled harmonic oscillators, described by 
the Hamiltonian Hl — J2 r e w H ( V r (. a H where r = 1,2 
labels the environments while I stands for the infinity 
set of oscillators whose frequencies are denoted by vD r g. 
ate (a>ri) represents the creation (annihilation) operator 



for the Ith mode of the rth environment. The action of 
these environments on the spin system is modelled by the 
interaction 



where we have defined the functions 



t+T 





"1 


2 lim 




T-S-0 


r 


2 lim 


"1 




r 



dx J dx' (ri(x)r r ,( x ')) 

1 r-t+r px 

- J dx J dx' (T r (x)Tl,{x')j 



(10) 



For the environments considered in this work it 
follows that: ^rt(i)rj,(t')\ = (T r (t)T r ,(t')) = 

0, (rt(t)IV(f)) = (nr)Srr> and (r r (f)rt,(f )) = 
((?i r ) + 1) 5 rr i, (n r ) being the thermal average excitation 
of the rth environment. These relations, of course, de- 
pend on the state of the environment. We observe that 
this master equation describes the Markovian evolution 
of a general time-dependent system, provided that the 
conditions i) and ii) of the last section are satisfied. 



IV. THE CHARACTERISTIC RELAXATION 
TIMES 

In this section, restricting us to the case of spin-1/2 
systems, we aim to derive the Bloch equations for the 
evolution of the magnetization components of N non- 
interacting spins. First, we obtain the Bloch equations 
from the Redfield formalism, relating the characteristic 
relaxation times with the properties of the associated 
classical stochastic environment. Next, computing the 
evolution of the average magnetization from the master 
equation formalism, we are able to link the characteris- 
tic relaxation times with the properties of the quantum 
environment. 



A. From the Redfield to the Bloch equations 



where T r (t) = '^2 l ^ri{t)ari and */ri(t) = 
j^ hr (t)ex.p[iA k (t)} with A k (t) being the phase 
factor coming from the transformation U' (t)Ola r tU(t) 
to the interaction picture. It is worth mentioning 
that the time-dependence of the system-environment 
coupling in the Schrodinger picture, r )^ hr {t) 1 comes 
from the assumption of a time-dependent spin system 
Hamiltonian, Hg(t). In fact, the coupling strength 
lri hr (t) leads to the decay rate of the master equation 
which plays the role of the time-dependent relaxation 
matrix in the Redfield equation ([6]). 

By inserting Eq. (10 1 into Eq. |9]) and performing 
the trace over the environments degrees of freedom, we 
obtain the master equation in the interaction picture 



da(t) 
dt 



Grr>(t)[Oi,a(t),O r }+H.c.}, 



(11) 



Let us consider here a spin-1/2 system placed in a con- 
stant magnetic field in z direction. The frequency gap 
between the two Zeeman levels defines the Larmor fre- 
quency lol (t) = lji (t) —loq (t), with uji (t) and ujq (t) rep- 
resenting the frequencies of the excited and the ground 
state, respectively. The modulation of these frequencies 
are due to some external influence, like an additional 
time-dependent magnetic field. The bare Hamiltonian of 
the spin-1/2 system is then given by H$(t) = ujl (*)Xz- 
The action of the environment over the system is mod- 
elled by the spin-lattice Hamiltonian 



H SL {t) = -l n J2\ q (t)I q , 



(12) 



where j n is the gyromagnetic factor, q labels the orthog- 
onal Cartesian directions {x, y, z}, X q (t) refers to the lat- 
tice stochastic fluctuation in q direction, and I q — u q /2 
stands for the spin (Pauli) operator. 
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For this system, the spectral density given in Eqs. ([8]) 
becomes 

Jkn,n'k'(ti^n'k') = E Iq ™ Iq k ® q (t,Sl n 'k>), (13) 

where I q = (k\ I q \n) and 

Q q (t, Sl n >k>) = I dt'e~^'^ T ° e in »'*' . (14) 



To derive Eq. ( 14 ) we have assumed isotropic stochastic 



fluctuations jS], by which 

\ q (t)\ q ,(t + t>) = VA'e-l 4 'l/ ro , 

being a mean value depending on the specific nature of 
the spin system and To the environment correlation time, 
measuring the rate of flips between the bath spins due 
to a specific anisotropic spin interaction (chemical shift, 
dipolar coupling, etc. ) [2D]. Note that we have assumed 
that the mean values of the coupling X q (t) are not affected 
due to the time-dependence of the system. This is quite 
reasonable since we are modelling the environment as a 
stochastic noise source. Remembering Sec. II, we have 



dr [oj k (t) - u n (t) 



which is just the integral of the Larmor frequency with 
a positive or negative signal, depending on the difference 
k — n (k,n = 0, 1). 

Next, by substituting Eq. (13) into Eq. we obtain 
the elements of the relaxation matrix 

iz*n,»'*' (*) = e { [ e ^' n "' k> y Ukn {t] 

1 

+e g (t,si kn )e in n'*'W 

3 

+S kn lfl^Q q (t, fi^Je^VW] } , (15) 

which enable us to compute the evolution of the mean 
value of the magnetization (Id) in an arbitrary d direc- 
tion: 



rkn rn ' k' 
1 8 



r dm 



By replacing Eqs. ^ and ( 15 ) into the r.h.s. of the last 
equation, we obtain 



d(I d 
dt 



= -iJ2l!T(m\[H s (t),E]\l) 

l,m 

+ E E 9 «(*> V ml )e-^T d nl (l\ [[I d , I q ] , E] 

q l,m 



Since, for spin-1/2 systems, |ri m /(i)| = Sl(£) (1 — 6 m i), the 
above equation can be separated into the longitudinal 
and transverse field components, 



dAjdj 
dt 



q—{x,y} l,m 

+ K s Tr{I z [[I d ,I,],H]}, 



where K q (t, Sl m i) — Q q (t, Sl m i)e^ lUml ^\ We stress that 
we have neglected the free-evolution term in the above 
equation because we are only interested in the effect of 
the environment induced dynamics. We also note that 
(m\ I x ( y ) \l) ^ if and only if m ^ I and (m| I z \l) ^ 
when m = I, which explains why the term K> x i y \(t, Sl m i) is 
a time-dependent function while k z is a constant. Writ- 
ing the last equation in terms of the longitudinal and 
transversal magnetizations, defined as M z — (I z ) and 
Mi = (I x ) x + (I x ) y, respectively, we obtain 



dM„(t) 
dt 

dM±(t) 



= -Re [ Kx (t, Si) + Ky(t, SI)} {M z (t) - M } , 



= --Re[/s x (t,n) + /s v (t,n) + 2/s,]M x (t), 

Mo = (I z ) eq = Tr{/ Z cr eg } being the equilibrium longitu- 
dinal magnetization. Now, comparing these results with 
the phenomenological Bloch equations 



dM z {t) 

dt 
dM ± (t) 
dt 



{M z (t)-M Q }. 



M ± (t) 

T 2 : 



(16a) 
(16b) 



the characteristic relaxation times, in terms of the time- 
dependent decay rates n q in the classical stochastic envi- 
ronment, is defined as 



Ti 
1 



= Re [n x (t,Sl) + K y (t,Sl)], 

= ^Rc [K x (t, SI) + K y (t, St)} + K z , 



(17a) 
(17b) 



Equations (16 1 and (17) show that, in contrast to the 



longitudinal rate Ti, the transverse decay rate T 2 is re- 
lated to an energy conserving process, affecting only the 
quantum coherence of the system. This fact justify the 
choice of the amplitude- and phase-damping channels for 
the quantum description of the spin system. It is worth 
mentioning that T\ can be controlled through the time- 
dependent parameter Sl(t) while T 2 can only be partially 
controlled since the decay rate k z does not depend on 
Sl(t). Finally, from Eqs. (17) we obtain the well known 
relation between both characteristic times 



1 



1 

2TY 



(18) 
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B. Prom the master equation to the Bloch 
equations 

In this section we consider the same system as before, 
but instead of a classical noise the spin system inter- 
acts with two quantum environments. In order to com- 
pute the evolution of the average magnetization from 
the master equation (11), which takes into account the 



amplitude- and phase-damping channels, we first address 
the decay rates T rr i (t) and Q rr i (t) defined in the end 
of Sec. III. When considering the amplitude-damping 
channel (r = a) we associate O a (0\) with the lowering 
(raising) spin operators whereas in the case of phase- 
damping (r = p) we define O p as the Hermitian number 
excitation operator. We then set O a — J_ and O p = I z . 
For both cases we set the mean values for the environ- 
ment operators (^ rl a r k) = (n r ,fe) Sik and, consequently, 

(a rk a\^j = ((n r . fe ) + l)6i k , (n r .k) being the thermal av- 
erage excitation of the fcth mode of the rth environment. 
Considering that the environment frequencies are very 
closely spaced to allow a continuum summation, such 
that J^i ~~ * (2 71 ") 1 dvJ r (v), J (v) being the spectral 
density of the environment, we obtain, for the amplitude 
damping case, the effective time-dependent decay rates 



F a (t) = ^e a (t), 



(19a) 
(19b) 



where, remembering the time-dependence of the 
system-environment coupling in the Schrodinger picture 



9„(t) = lim 



T-fO 



t + T pX fOO 

dx I dx' 



dvJ a (v) 



X <yf chr (u, x')^ chr {u, x ) e iH*')-n(x)+u(x-x')] 



To obtain the effective decay rates in Eq. ( 19 ) it was as- 



sumed, as usual, that the thermal average excitation of 
the environment modes vary slowly around the range of 
variation of the spin system frequency. This is a good ap- 
proximation when the environment is in a thermal state, 
as the present case [9]. 

For the phase-damping channel, the effective decay 
rates are given by 



(20a) 
(20b) 



where we defined 



Q p (t) = lim 



t+T 



dx j dx' 



duJ p (v) 



x ^ chr (iy,x')^ chr (^x)e^ x '- x ^ 



Due to the diagonal system-environment coupling asso- 
ciated with the phase-damping case, we are able to com- 
pute Q p (t) without having the explicit form of fl(t). In 
fact, assuming that the spectral density J p {v) as well as 
the system-environment coupling J p (v, x') vary slowly 
around v = 0, wc obtain the time- independent parameter 



0„ — lim 

y T->0 



4.(0) 

T 



t + T 



[^ Chr (0,x)} dx 



- J vU P \ ■ 

From this result we see that, in contrast to the case of the 
amplitude-damping channel, the decay rates T p and Q p 
do not acquire a time dependence due to the modulation 
of the system frequency, resembling the result obtained 
in the Redfield formalism. Finally the master equation 



( 11 1 becomes 
da(t) (i 



dt 



2tt 



-0 o (t) [J-<r(i),/-+ 



(("») + 1) 
2tt 

|(2<n p ) 



o (t) [I+tr(t)J-] 

+ l)[IMt),I z ]+H.c. (21) 



Computing the evolution of the mean value of the mag- 
netization (1^) in an arbitrary d direction we obtain for 
the longitudinal and transversal magnetizations 



dM z Rc9 a (t) 

dt 7T 

dM ± _ [Ree Q (t) 
dt ~ [ 2^ 

Re6 ,„ , , 



(2 (n«) + l)Af„ 
(2(n„) + l) 



1) 



Mi 



As in the preceding subsection, we compose these equa- 
tions with the Bloch equations (16), to obtain 



1 _ Re0 a (t) 



1). 



T\ it 
1 RcQ a {t) 
T 2 ~ 2n 
and, consequently, the relation 
1 1 Re9 p 
% ~ 2J\ + 7T 



(2 (n,) 

(2 (n a ) + 1) + (2 (n p ) + 1) ; 



RcQ r 



(22a) 
(22b) 



(2 (n p ) + 1) . 



which has the same structure as Eq. (18), obtained by 
the Redfield formalism. From Eqs. (17) and (22) we can 
make the identifications 



Re6„(i) 



7T 

Re6 



(2 (n a ) 



1 



Re[n x (t,ty 



E (2<n p ) 



1) = K, 



■Ky(t,n)], (23a) 
(23b) 



These equations show the connections between the semi- 
classical and quantum approaches to open system dy- 
namics. In the next two sections we will construct the 
Kraus and the phenomenological operators for the time- 
dependent system studied in this section. 
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V. THE OPERATOR-SUM REPRESENTATION 

It is well known that every transformation that is given 
by a completely positive map admits a representation of 
the form ffl] 



(24) 



with the Kraus operators Ek (t) satisfying the following 
relation 



Y,E$(t)E k (t) = l. 



Our goal in this section is to construct the operators 
Ek (t) for both channels studied in the last section. To 
achieve this we will consider the density operator evo- 
lution equations, which follows from the Redfield or the 
master equation formalisms. Then, we compare these 
equations with those shown in Eq. ( 24 1 to obtain the 



time dependence of the Kraus operators. 

In the next two subsections we will adopt the basis de- 
fined in Sec. IV, that diagonalize the spin operator I z , 
{|0) , |1)}. Since both channels studied here are indepen- 
dent, let us then consider each one of them separately. 



A. Phase-damping channel 

Considering only the phase-damping (6 a = 0), the 
master equation (21 1 leads us to the following set of dif- 



ferential equations satisfied by the elements of the density 
operator 



daujt) 

dt 
duoo(t) 

dt 
d<r w (t) 

dt 
dapi (t) 

dt 



-r p o-io(t), 
-r p cr i(t), 



(25a) 
(25b) 
(25c) 
(25d) 



where we have defined T p — 2Re \T P + G P ] — 
ReQp (2 (tip) + 1) /tt. Note that, as expected, the popu- 
lations are not affected by this noisy channel. We assume 
that the Kraus operators for this case are given by [5] 



1 

-1 



(26a) 
(26b) 



with p (t) being the parameter to be determined. Starting 
from the initial density operator 



a(0) = 



o 



'01 



'00 



and substituting Eq. (26) into Eq. (24 1, we thus obtain 

(28) 



fau[t) <rio(f) 
V?oi(<) °oo(t) 



'ii 



l-2p)a° 



(l-2p)< 



01 



'00 



By imposing that the time derivative of the elements of 
the above evolved density operator must be identical to 
those in Eqs. (25), we derive the following differential 



equation for the parameter p (t) 
dp r 



dt 



?(2P-1), 



with the initial condition p (0) = 0, arising from the fact 
that the Kraus operators must reduce to the identity at 
the initial time. The solution for p (t) is thus given by 



v (*)= *{i-exp[-ry]}, 



(29) 



which finally defines the Kraus operators in Eqs. (126 



B. Amplitude-damping channel 

For the case of amplitude damping we assume that the 
Kraus operators are given by [5] 



El 



El 



El 



Ht 







y/l-a(t) 
y/atfj 



,0 

V 1 ~ 7t 







^l-a(i) 
1 



^^{Jk) o)' 



(30a) 
(30b) 
(30c) 
(30d) 



where jt = exp[—(3E]/Z is the Boltzmann factor, E 
being the energy gap of the spin- 1/2 levels and Z = 1 + 
exp [—/3E] the partition function. By analogy with the 
preceding subsection our aim is to obtain the differential 
equation obeyed by the parameter a (t). To this end, we 
compute the time evolution given by Eq. ( 24 ) , to obtain 



o-n(i) <Tio(t)\ _ fj T a + q (l - ag ) 



(1 - 7t) a + qa^ 



with q = 1 — a. By imposing that the differential equa- 
tions derived from the above density operator must be 
identical to those obtained from the master equation 
(21 ), with <d p = 0, given by 

dan(t) 



dt 
daiojt) 

dt 
da i jt) 

dt 
d<J 00 (t) 

dt 



2Re { T a (t)<7 00 (t) - Q a (*)} , 
-{K(t)+0a(t)}a 10 (t), 
-{F a {t)+g* a {t)}a 01 {t), 
-2Re {T a (t)a QQ (t) - G a (t)a n (t)} . 
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we finally obtain 

[lT ~ <&] ^ = - 2Re{T a (t) + g a (t)} ( lT a + qa ^) 



2Re{F a (t)}. 



da 

lit 
da 

Itt 



a° w =2{T a (t)+g:(t)}qa° 10l 
a° Q1 =2{F* a (t) + g a (t)}q<j Q1 . 



The fourth equation being identical to the first one. 
Adding the last two equations and remembering that the 
initial conditions for the elements of tr(0) are arbitrary, 
the following set of differential equations must be satis- 
fied 



7t da 



da 
~di 



=Bje{^ a (t)+g a (t)}( 7T a + qa^) (31) 
+ Re{J- a (i)}, 



--2Re {T a {t)+g a (t)}q, 



(32) 



together with the condition a(0) = 0. The solution of 
Eq. ( 32 ) is readily obtained as 



a (t) = 1 — exp 



-2 f Re{T a (r)+g a (T)}dT 
Jo 



(33) 



and its substitution into Eq. (31 1 leads to the relation 
( 7t - 1) Re { F a {t)} + 7T Re {g a (t)} = 0, 



which, together with Eqs. (19), results in the identity 

/ \ 7t 1 



1 - 27 T 



°J E - 1 ' 



Equation ([33]) satisfies both Eqs. (31) and (32), thus 
being the desired solution. 

It is now straightforward to see the connection between 
both ab-inito approaches, the Redfield and the master 
equation, with the operator sum representation. Back 
to the identifications in Eqs. (23) and the definition of 
the relaxation times in Eqs. (22), we can rewrite the 



time dependence of the Kraus operators for the phase- 
damping channel as 



Pit) 



exp 



K z t 

~2~ 



and for the amplitude-damping channel as 



- (<) = 1 — exp 



1 



o Zi(r) 



It is interesting to point out that within the operator 
sum formalism the distinction between the decay of the 
coherences and the T 2 decay in NMR systems becomes 
evident. While the density matrix coherences decay 
is completely independent from the T\ relaxation time, 



there is an intrinsic dependence of T 2 with T\ through 
Eq. (18). Thus, one can conclude that the phenomeno- 



logical T 2 does not reflect only the decay of the quantum 
coherences of the system. 

Before applying this formalism to a specific situation, 
let us make some remarks regarding a very closely related 
technique, the phenomenological-operator approach [TT] 
also introduced to simplify the treatment of dissipative 
quantum systems. 



VI. THE PHENOMENOLOGICAL-OPERATOR 
APPROACH 



The phenomenological-operator approach is a tech- 
nique equivalent to the operator sum representation, but 
taking explicitly into account the state of the environ- 
ment together with those of the open quantum system. 
From the phenomenological operators we automatically 
derive the Kraus operators and vice-versa. 



A. Phase-damping channel 

First, let us consider the coupling of the spin- 1/2 states 
to a surrounding phase-damping environment, which can 
be described by the map 



|0)|O^|0)f ( p) |£> 
|1)|£)->|1)7#>|£) 



\Vt$\e) 



(34a) 
(34b) 



where \£) denotes the initial state of the environment and 
the operators 7", acting on this state, account for the 
system-environment coupling. Since a phase-damping 
channel does not exchange energy with the system, we 



obviously have the identity operator 77, 



(/') 



Re- 



garding the excited initial state |1), it will remain as 
such, with or without an additional phase shift relative 
to the ground state |0), due to the action of the envi- 
ronment. In the later case, we must impose (after the 
computed master equation (21) with 6 a = 0), the decay- 



ing probability T^f = e _rp *l, with the above defined 
rate T p — 2Re [F p + g p ] , remembering that the time- 
dependent frequency does not lead to a time-dependent 
relation rate T p . In the former case we have T 10 



(t) ^a] + OjJ , with fj(t) giving the probability am- 
plitude for environment, described by the creation and 
annihilation operators and aj , respectively, to induce a 
phase shift on the excited state of the system. Assuming 
that the state of the environment is modified such that 
(£|7^|£) = 0, we obtain, after normalization of the 
state vector |1)|£), the relation £\ |fj(*)| 2 = 1 - e" 2 ^*. 



It is straightforward to verify that the map in Eq. (34) 



leads exactly to the density operator ( 28 ) derived for the 
phase-damping process. 
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B. Amplitude-damping channel 



VII. COHERENCE CONTROL 



Before addressing the non-zero temperature case, we 
first consider the coupling of the spin-1/2 states to an 
amplitude-damping environment at T — K, described 
by the map 



|0)|£> 



|o>r 



o ( o) |£>, 



\l)t^\£) 



\0)T$>\£ 



With the environment in the vacuum state we obviously 
have the identity operator Tqq^ = 1, and after the com- 
puted master equation ( |21[ ) with O p = 0, we must impose 

that 7#° = exp [-2/ t Re{J- a (r) + G a (T)}dr 
time-independent system we obtain the usual solution 



1. For a 



T-O) 
'11 



e - r -'l, where T a = 2Re [T a + G a ]. For the op- 



erator 7^ associated with the excitation of one of the 
infinite environment modes, we have T~io = Sj(^) a ]) 
with gj(t) giving the probability amplitude for the ex- 
citation of the jth oscillator mode of the environment. 
After normalization of the wave vector we obtain 

£il(&(*)l a = l-e- ar -*. 

By turning our attention to the case of a non-zero tem- 
perature environment, we can write the extended map 



|0>f o ( o a) |£) 
\l)t£ ] \£) 



\0)T^\£), 



where, now, instead of the identity operator — 1, we 
must account for the probability of excitation of the sys- 
tem due to the environment background photons. Since 
the operator 7~qq^ is associated with an event at which the 
environment is not excited, it must remain proportional 
to the identity. Moreover, as far as the probability for 
the system to remain in the ground state must decrease 
in a rate proportional to the environment temperature, 



we naturally impose that Tq 



(a) 



00 



)7 T 1. 

From the above assumption for Tqq^ we straightforwardly 
obtain for 7q " = £ • hj(t)a,j, the relation J2j |hj(i)| 2 — 



,-2r„( 



)~/Ti with gj(t) giving the probability am- 



plitude for the system to be excited by the jth oscillator 
mode of the environment. Regarding the operator 7ii, 
we know that it must be also proportional to the identity 
since the environment must remain unaffected. However, 
for the case of non-zero temperature, the probability for 
the system to remain in the excited state must decrease in 
a rate smaller than the exponential decay factor e _2ra */ 2 
coming from an environment at absolute zero. Moreover, 
the equilibrium probability must depend on the ther- 
mal average photon number (n a ), such that we impose 



-2T a tt 



(1 



-2T a tt 



the complementary operator 7^ = Ylj "j 



Consequently, for 



(t)a]j we ob- 



tain E ; - Mi)| 2 = y /(l-e-r t )(l-rr) 1 h,-(t) being the 
probability amplitude for the excitation of the jth oscil- 
lator mode of the environment at non-zero temperature. 



In Ref. 1 1 6 j it was proposed a method to circumvent 
the decoherence process of a non-stationary system un- 
der amplitude-damping channel. In that work, the au- 
thors focused on the state protection of a cavity mode 
whose modulation of the frequency ui(t) was engineered 
through the atom-field interaction. The master equation 
approach was used to investigate the dynamics of the cav- 
ity mode, assuming its interaction with the environment 
to be proportional to 



(36) 



[u{t)-»Y+e 



v being the continuous frequency of the environment and 
the parameter £ accounting for the Lorentzian sharp- 
ness of the coupling around the frequency u)(t). Thus 
the Lorentzian coupling, which is justified in the system- 
environment weak-coupling regime, "follows" the evolu- 
tion of the frequency of the system, as expected under 
the sudden coupling approximation. 

The modulation of the frequency was engineered to be 
of the form 



uj(t) = loq + xsin^i, 



(37) 



with wo being the static frequency of the cavity mode. 
The condition C/ w o ^ 1> easily achieved within typical 
experimental conditions as in NMR and cavity quantum 
electrodynamics, define the adiabatic modulation of the 
frequency. Out of this regime, when £ > ujq we reach 
the regime of the Casimir-like effect [22], where the de- 
coherence mechanism of the cavity mode are completely 
distinct from the one discussed here. 

Within the condition C/ w o < 1, wc show that the con- 
trol of the decoherence process is achieved by means of 
the two parameters 



n 



e = 



C ' 

£ r 



X 



X 



(38) 
(39) 



where Tq is the natural decay rate of the cavity mode and 
we have assumed, as it is expected, that £ ~ Tq. It was 
demonstrated in Ref. |16j that a significant attenuation 
of the decoherence occurs when both of these parame- 
ters are smaller than unit. This is seen from the derived 
time-dependent decay rate of the cavity mode, which gets 
weaker as one or both parameters r\ and e decreases. The 
physical reason for this can be seen as follows: The char- 
acteristic time interval for an appreciable action of the 
environment over the stationary system is around Tq 1 . 
However, when the frequency of the system changes con- 
tinuously, its rate of variation (proportional to £) plays a 
crucial role in the effective coupling between the system 
and the environment. Remembering that this coupling 
occurs around ui(t), in a region defined by the Lorentzian 
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sharpness £ [see Eq. ( 36 )] , a rate of variation C signifi- sion for D 



cantly larger than To, such that 77 1 , makes difficult 
an effective action of the environment over the system 
since their interaction time is reduced proportionally to 
rj. Otherwise, when Q is smaller than To, an effective 
action of the environment takes place, inducing the re- 
laxation of the system before a significant change of its 
frequency. By its turn, the role of the amplitude of the 
oscillation x is to trigger the action of the rate of vari- 
ation £. In fact, when the amplitude x 1S smaller than 
the Lorentzian sharpness £, the non-stationary system 
does not leave the region (in frequency space) of its ef- 
fective coupling with the environment, thus decaying as 
a stationary system, whatever the value of £. However, 
when x is larger than £, the effective system-environment 
coupling moves to different regions of the spectrum, thus 
triggering the action of the rate of variation £ as described 
above [16] . 

Now, we apply the same idea to the case of a spin-1/2 
system considering the Redfield formalism to treat the 
decoherence process. The same conclusions of Ref. [TS] 
are obtained here, but without imposing any specific form 
for the system- environment coupling, as in Eq. (36), since 



a Lorentzian time-dependent effective decay rate for the 
spin-1/2 system automatically appears from the Redfield 
formalism applied to NMR systems. As far as we do 
not have to define the function in Eq. (36), we stress 



that the parameter e, which takes place explicitly in the 
effective time-dependent decay rate of the cavity mode in 
Ref. [H], does not appear in the present spin-1/2 case. 
Instead, we must use the ratio x/C also defined in Ref. 
[16j . weighting the contributions of parameters 77 and e. 

The spin-1/2 system, interacting with the amplitude- 
damping environment, is described by the Hamiltonian 
Hs(t) — ljl {t) I Z i with u>l (t) being modulated as in Eq. 
(37). We note that such a frequency modulation beras 



no conexion with To see how the protocol works, con- 
sider the solution of the differential equation governing 
the evolution of the longitudinal magnetization [see Eq. 
(16a)], which decays with the effective rate given by 



T\ (t) defined in Eq. (17a). Therefore, to circumvent de- 
coherence we have to make T\ (t) grater than its static 
value (equivalent to T^ 1 in the previous case of a 
damped cavity mode). This fact can also be directly 
seen from the operators in Eqs. (30), which reduces to 
the identity when a{t) — >• 0, i.e., when D (t) — » 0. By 
considering an isotropic and homogeneous environment, 
Eqs. (14), (17a), and p7|) leads to the following expres- 



r T2 

D (r) -2 ( 7 nAT°) 2 J dr 2 J dn exp |-^| 


x cos jlfWLTi + 2^ sin 2 f j , (40) 



where we have performed the ch ange of variable r = t /T® 



and r\ = (T^C) as in Eq. (38). The gyromagnetic fac 



tor 7„ and the lattice stochastic fluctuation in q direction 



A = X„ is defined in Eq. ( 12 ) 



In Fig. 1 we plot the function a(r) = 1 — c^ D ^ against 
the dimensionless time r, for some values of 77 and the 
ratio x/C- As can be seen, the relaxation decay rate for 
the stationary case x = 0, the black solid curve, is sig- 
nificantly attenuated when the control parameter 77 de- 
creases and/or the ratio x/C increases. The red-dashed 
and the blue-dotted curves correspond to the pairs of 
values (ij — 0.1, x/C — 1) and (v = 0.1, x/C = 10), re- 
spectively. The curves corresponding to the pairs (77 = 
0.01, x/C — 1) an d (77 = 0.01, x/C = 10) are very close to 
the red dashed and blue dotted lines, respectively, and, 
therefore, they are not shown in the figure. From this 
fact we conclude that the role of parameter x (the am- 
plitude of the modulation) in the decoherence control is 
more effective than the role played by C (the modulation 
frequency). This result is in perfect agreement with the 
one obtained in Ref. |16j . An important observation is 
that this protocol is completely independent of the tem- 
perature of the environment, as well as of the initial state 
of the system. 

A possible implementation of the proposed scheme 
could be realized in an NMR experiment as discussed 
below. In addition to the static Zeeman field B which 
defined the Larmor frequency u>l = 7n-Bo, another paral- 
lel time-dependent field B(t) is applied to modulate the 
frequency of the spin-1/2 system. This additional field 
may be furnished by a Hclmholtz like coil surrounding 
the probe, traversed by a tailored current which provides 
the time-dependent component j n B(t) of the frequency. 
By imposing the time-dependence of the auxiliary field 
to be of the form 

ln B(t) = X sm (Ct), 



we obtain the frequency modulation given by Eq. (37) 



To circumvent decoherence we must have 77 <§C 1 [see Eq. 
(38) and Fig. 1], which implies that the modulation fre- 



quency must obeys 



c»(7fT 



(41) 



This relation is in agreement with the adiabatic condi- 
tion, since we have, in general, T% 3> w^ 1 . The other 
condition for the protocol to work is that the amplitude x 
of the oscillation of the system frequency must be greater 
than the spectral sharpness of the coupling which, as said 
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FIG. 1. (Col or o nline) Plot of the decay function a(t) de- 
fined in Eq. ( 33 1 for some values of the control parameters 



r\ and x/C The black solid line represents the uncontrolled 
case, where \ = 0. The other two shows the effective control 
of the relaxation process due to the modulation of the sys- 
tem frequency. The vertical green dash-dotted line marks the 
relaxation time for the static (uncontrolled) case. 



before, is expect to be of the order of (T-f) , such that 
X^iny 1 . (42) 



Besides the conditions (HT1) and (42 1, the magnetic field 



B{t) must also obeys the adiabatic condition [see equa- 
tion bellow Eq. 0)], which leads to \B(t) \ < \B \. This 
conditions can be easily attained in a typical NMR ex- 
periment. Considering a Hydrogen nuclei at B ~ 10 
T (uj l ~ 400 MHz) and T x ~ 1 s, we obtain, together 
the following range for the magni- 
field B(t) for the protocol to work, 



with condition (41 



tude of the contro 
i.e., (jnT^y 1 « \B(t)\ « \B \: 



lQ- b T < \B(t)\ < 10 T. 

Note that the NMR setup is appropriate for this kind 
of experiment due to the fact that the relaxation time is 
relatively large compared with other platforms. This fact 
permit us to realize the experiment in a high-temperature 
environment. It is also important to observe that despite 
in the NMR experiments both amplitude- and phase- 
damping are present, there exist experiments where only 
the amplitude-damping are probed, so they can be used 
to check out our proposals. Therefore, we can only con- 
trol the relaxation time T\, but not Ti. 



VIII. FINAL DISCUSSIONS 

We have presented here a unified view of the semiclas- 
sical Redfield formalism and the quantum master equa- 



tion approach for a time-dependent spin system. Focus- 
ing on a spin-1/2 system, we shown the equivalence be- 
tween both approaches through the fact that they lead to 
the same Bloch equations and, consequently, to the same 
characteristic longitudinal T\ and transversal T2 relax- 
ations times. We verified that only T\ is affected by the 
time-dependency of the system frequency and built the 
Kraus and the phenomenological operators for the spin- 
1/2 system under both the amplitude- and the phase- 
damping channels assumed within the master equation 
approach. 

As an application, we revisited a protocol to circum- 
vent relaxation and, consequently, the coherence control 
of a non-stationary system 16J. In contrast to the pro- 
tocol in Ref. [16] . in the present case, we do not imposed 
a functional form for the system-environment coupling, 
which emerges naturally from the Redfield semiclassical 
formalism. The coherence control of the spin-1/2 sys- 
tem was demonstrated by enlarging the longitudinal re- 
laxation time through the modulation of the system fre- 
quency. We stress that the protocol in Ref. [16], trans- 
lated here for the spin-1/2 system is different from the dy- 
namical decoupling methods presented in literature [15] . 
where one must interfere in the system in time scales less 
than the bath correlation time. The frequency modula- 
tion technique takes advantage of the pre existing natural 
frequency of the system, adding to it a small amplitude 
to achieve such time scales. We also discussed the im- 
plementation of this protocol to control the longitudinal 
relaxation time in the NMR context. 

It is worth stressing that the development present 
here applies for Markovian environments and adiabatic 
modulation of the required time-dependent frequency. 
Therefore, an extension of the present development for 
non-Markovian environments as well as for interacting 
time-dependent systems would be desirable, since non- 
Markovian environments is present in many promising 
platforms for quantum information processing such as 
photonic crystals [23] ■ Moreover, time-dependent inter- 
acting system would considerably enlarge the perspective 
of the present work for quantum information purposes. 

Although the entanglement is not usually present in 
the NMR system [53] , there exists classical and quantum 
correlations that could be useful for some quantum infor- 
mation protocols [25]. The decoherence control exposed 
here can then be directly applied to protect these correla- 
tions from the action of the amplitude-damping channel. 
Finally, we stress that the theory presented here is com- 
pletely general and can be applied to different platforms 
for quantum information processing such as cavity quan- 
tum electrodynamics, trapped ions, quantum dots, etc. 
However, for this statement to be valid we must assume 
that all the rotations needed to implement quantum in- 
formation processing occur at the beginning or at the end 
of the experiment; in between them, the system evolves 
only under the action of the environments. 
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